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to the recent full-QCD lattice data. We find good agreement for the case of the Spectral Quark 
Model, which builds in the vector-meson dominance for the charge form factor. We derive a simple 
relation between the gravitational and electromagnetic form factors, holding in the considered quark 
models in the chiral limit. The relation implies that the gravitational mean squared radius is 
half the electromagnetic one. We also analyze higher-order quark generalized form factors of the 
pion, related to higher moments in the symmetric Bjorken X-variable of the generalized parton 
distribution functions, and discuss their perturbative QCD evolution, needed to relate the quark- 
model predictions to the lattice data. The values of the higher-order quark form factors at t = 0, 
computed on the lattice, also agree with our quark model results within the statistical and method 
uncertainties. 
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I. INTRODUCTION 

Form factors carry basic information on the extended 
structure of hadrons as they correspond to matrix el- 
ements of conserved currents between hadron states. 
There exist abundant experimental data concerning the 
charge form factor, related to the electromagnetic current 
and providing the distribution of charge in a hadron. The 
gravitational form factors, related to matrix elements 
of the energy-momentum tensor [l[ in a hadronic state 
and thus providing the distribution of matter within the 
hadron, are not experimentally known. This is because 
whereas the electromagnetic interactions are probed by 
structureless electrons and mediated by one-photon ex- 
change, the parallel one-graviton exchange is extremely 
weak and impossible to measure. Recently, however, 
these objects were determined ab initio and with suf- 
ficient accuracy in full-QCD lattice simulations by the 
QCDSF/UKQCD collaboration [2, 3]. In this paper we 
confront the predictions of chiral quark models for the 
pion gravitational form factors and higher-order gen- 
eralized form factors to the lattice determination of 
Refs. U, 3]. We find that in the Spectral Quark Model 
(SQM) which is a variant of a chiral quark model with 
built-in vector meson dominance, the agreement with the 
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lattice data is especially good. In particular, a slower fall- 
off with the large space-like momenta than for the case 
of the electromagnetic form factor is found. We also per- 
form the calculations in the Nambu-Jona-Lasinio (NJL) 
model with the Pauli-Villars (PV) regularization, where 
the agreement is not as good as in SQM. In addition, 
we find that the higher-order generalized form factors 
(see below for a definition) at t = provided by the 
full-QCD lattice simulations are properly reproduced in 
chiral quark models. 

The form factors are related via sum rules to more gen- 
eral objects, the generalized parton distributions (GPDs) 
of the pion (for extensive reviews see e.g. H H, @, H H 
Hoi. IT1II12I Il3j] and references therein). Experimentally, 
the GPDs of the pion are elusive quantities, as they ap- 
pear in exclusive processes which are difficult to measure, 
such as the deeply virtual Compton scattering or the hard 
electro-production of mesons. Recently, it has been sug- 
gested to study instead the deeply virtual Compton scat- 
tering on a virtual pion that is emitted by a proton [l4| 
under the operating conditions which will first be met 
after the energy upgrade at TJLAB. This will eventually 
set important constraints on the pion GPDs. 

Despite their fundamental and general character, 
GPDs are genuinely defined in the Minkowski space, thus 
hindering direct determinations on Euclidean lattices. 
The moments of GPDs in the X variable (see Sect. [TTJ) 
form polynomials in the £ variable, with coefficients de- 
pending on the t- variable only. The lowest moments yield 
the standard electromagnetic and the gravitational form 
factors, while higher moments are known as generalized 
form factors. These are useful quantities which may be 
computed directly on the lattice. Presently, apart from 
the values at t = 0, the generalized form factors are 
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not known even from lattice simulations due to insuf- 
ficient statistics. In this work we make predictions for 
the first few generalized form factors. Since these ob- 
jects do not correspond to conserved currents, they evolve 
with the QCD scale as they carry anomalous dimensions. 
We undertake this perturbative analysis at leading order 
and show how the generalized form factors evolve with 
the scale. We use the techniques described in detail in 
Ref. [3. 

Chiral quark models have proved to properly describe 
the essential features of the pionic GPDs and related 
quantities. The special case of the parton distribution 
function (PDF) has been analyzed in the NJL model 
in Refs. [IE [lj, [l8[ . The diagonal GPD in the impact 
parameter space was obtained in [19| . Other calcula- 
tions of the pionic GPDs and PDFs were performed in 
instanton-inspired chiral quark models [IE [U [22|, [23|, 
US HE El- In the effective quark models the GPDs 
have been analyzed in [IS HH, 0, [2E [3(| • Studies pay- 
ing particular attention to polynomiality were carried 
out i n [2lL [3ll I32I ] , which proceeded via double distribu- 
tions [33j ] . The same technique was applied in Ref. [H[, 
which allowed for a simple proof of polynomiality in chi- 
ral quark models. Also, other formal features, such as 
support, normalization, crossing, and the positivity con- 
straints [lEUE], are all satisfied in the chiral quark model 
calculation of Ref. [HI]. Moreover, the obtained analytic 
expressions have a rather non-trivial form which is not 
factorizable in the t-variable. The parton distribution 
amplitude (PDA), related to the GPD via a low-energy 
theorem l3ol was evaluated in Refs. [13, HE HE HE [41] 
(see Ref. [42| for a brief review of analyses of PDA). 
The authors of Ref. [3] provide standard electromagnetic 
and gravitational form factors based on a phenomenolog- 
ical model factorization assumption which is generally 
not satisfied by our dynamical field theoretical calcula- 
tion [TEl ]. Finally, the related quantity, the pion-photon 
transition distribution amplitude (TDA) H Q gl [i| 
has been obtained in Refs. [13, H HE HE1I1 • 

In the present work we find that in the considered class 
of models a simple relation between the gravitational and 
charge form factors holds (see the Appendix). The essen- 
tial element for the proof is the existence of the spectral 
representation of the quark model, which is the case both 
for SQM and NJL with the Pauli-Villars regularization. 
The relation implies that in the considered models and 
in the chiral limit the gravitational mean squared radius 
is half the electromagnetic one. 

In Ref. [l5| we have stressed the relevance of the QCD 
evolution for the phenomenological success of the chiral 
quark models in the description of the experimental and 
lattice data for the PDF and PDA of the pion. Indeed, 
the evolved results for the valence PDF compare very well 
to the Drell-Yan data from the E615 experiment [52| at 
the scale of 4 GeV, and to the transverse lattice results 
[5^ at lower scales, ^0.5 GeV. Similarly, for the case of 
the PDA the QCD evolution leads to a fair description 
of the E791 dijet data [1J] at the scale of 2 GeV, and of 



the transverse lattice data [53[ at the lattice scale. The 
comparison is presented in Figs. 8-11 of Ref. [l5[. 

The above-mentioned success of the chiral quark mod- 
els in describing properties related to the pionic GPD 
allows us to hope that also their other aspects, such as 
the form factors, including the generalized ones, can be 
reliably estimated in these models. Certainly, the elec- 
tromagnetic form factor, as one of the most basic quan- 
tities, has been promptly computed in all chiral quark 
models of the pion. Note that these evaluations assume 
the large- N c limit, thus the effects of the pion loops, im- 
portant at low momenta, are absent. We note that the 
NJL model leads to somewhat too small electromagnetic 
radius [41], while in SQM, which incorporates the vec- 
tor meson dominance, one may fit the data accurately. 
Some predictions of the NJL model for the generalized 
form factors of the pion have been presented by one of 
us in Ref. (55j , however, the data away from the physical 
pion mass have been used in that work. 

The outline of the paper is as follows: In Sec. [11] we re- 
view the necessary formalism providing basic definitions 
and notation. The essential information on the two chi- 
ral quark models used in our work is provided in Sec. IIII1 
while Sec. IIVI shows their predictions for the electromag- 
netic and gravitational form factors. In Sec. [Vj contain- 
ing our main results, we compare our predictions to the 
full-QCD lattice data from Refs. Then in Sec. EJ 

we pass to the results for the higher-order form factors, 
where the QCD evolution effects are incorporated. 



II. BASIC DEFINITIONS 

In this Section we review the basic concepts necessary 
for our analysis and introduce the notation. 

Throughout this paper we work for simplicity in the 
strict chiral limit, 

m„ = 0. (1) 

Since the lattice data of Refs. @, 0] perform the extrap- 
olation to the physical pion mass which is small, the as- 
sumption (fl]) is appropriate. Nevertheless, the exten- 
sion to the physical pion mass in chiral quark models 
is straightforward. In any case, it is worth mentioning 
that although an attempt was made [E H[ to incorpo- 
rate chiral logarithms from Chiral Perturbation Theory 
(%PT) to the one-loop order [HE [57l.l58l (for a review see 
e.g. Ref. [59j]), as described in Ref. |60l |, the data did not 
exhibit their presence when all other uncertainties were 
considered, suggesting instead a linear extrapolation. 

In this paper we will deal with standard gravitational 
and vector form factors. The corresponding quark oper- 
ators are 

9^= J2 9(^(7^ +7^)^), (2) 

q=u,d,.. . 
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and 



E 9»y7^)- 

q—u,d,.. . 



(3) 



respectively. The gravitational quark form factors of the 
pion . Oi and O2, are defined through the matrix 
element of the energy-momentum tensor in the one-pion 
state, 



(ttV) I e""(o) 1 w a ( P )) = 

h<* ^u q 2 _ q ^ Ql (q 2 ) + 4W8 2 (<7 2 )] , 



(4) 



where P = \{p' + p), q — p' — p, and a, b are the isospin 
indices. The gravitational form factors satisfy the low- 
energy theorem 9i(0) - 6 2 (0) = 0(m\) [gl|. In the low 
momentum and pion-mass limit one can establish contact 
with xPT in the presence of gravity [6l| . 



e 1 ( g 2 ) = i + ^(4L 11 



e 2 (<z 2 



1 



P 
~P 



To - (ill — ^13) 



L\2 



(5) 



where / = 86 MeV is the pion weak decay constant in 
the chiral limit and in. 12, 13 are the corresponding grav- 
itational low energy constants. The calculation of the 
effective energy-momentum tensor as well as the low en- 
ergy constants for chiral quark models has been carried 
out in Refs. [62, 63]. The vector form factor is defined as 

(7rV)l^ b (0)k c (p)) = e abc {p^+p^)F v {q 2 ), (6) 
At low momentum one has 



F V (q 2 ) = 1 + 



2q 2 L 9 
P 



(7) 



where Lg is a low energy constant of xPT [58(] . 

The generalized quark form factors of the pion (here we 
take 7r + for definiteness) are defined as matrix elements 
of quark bilinears accompanied by additional derivative 
operators, 



iD^- l} u(0) \tt + ( P )) 



(8) 



■ r 



[ P W ...P^ } 2- k A nk (t), 



k=2 



where A n k(t) are the generalized form factors with n = 

1,2,... and k = 0, 2, . . . , n. The symbol D is the QCD 

covariant derivative, D = h(D — D), and {...} de- 
notes the symmetrization of indices and the subtraction 
of traces for each pair of indices. The factor of 2~ fc is 
conventional and makes our definition different than in 
Refs. [2j, |3J]. For n = 1 and 2 we have 



A 10 (t) = F v (t), 

Aao(t) = ie?(i), A 22 {t) 



(9) 



where Fy is the pion electromagnetic form factor, while 
Ql denote the quark parts of the gravitational form fac- 
tors of Eq. flU). 

Chiral quark models work at the scale where the only 
explicit degrees of freedom are quarks, while the gluons 
are integrated out. Thus, the gluon form factors of the 
pion vanish, 



Ank M — (quark — model scale) . 



(10) 

When evolution to higher scales is carried out [l5l |. non- 
zero gluonic moments are generated. Note that the above 
mentioned soft-pion theorem &i(0) — 62(0) + 0{m 2 T ) 
applies to the full trace of the energy momentum tensor, 
and consequently does not apply to their quark or gluonic 
contributions separately. The condition (|10[) translates 
then into ^20(0) = —^22(0) in the chiral limit at the 
quark model scale. 

Now we pass to the definition of the GPDs, whose 
moments are related to the form factors listed above. 
The used conventions are the same as in Ref . [HI , where 
also all the details of the calculation of the GPDs can be 
found. The kinematics of the GPDs can be read off from 
Fig. [T] The adopted notation is 



0, q> 



Q 



0, p ■ n = 1, q ■ n 



-2p-q = t, 

<, (11) 



with the null vector n defining the light cone. The two 
isospin projections of the quark GPDs of the pion, isos- 
inglet (singlet) and isovector (non-singlet), are defined 
through the matrix elements of quark bilinears, with the 
quark fields displaced along the light cone, 



Sab'H 



1=0 



(•<■■<;■/)= / ^ ixp+z ~ 



(12) 



x (n b (p + g)|^(0) 7 • n^(z)\n a (p))\ z+=0 z±=Q , 



i^abT-L 1 1 (x,(,t) 



dz 

47T 



(13) 



x (ir b (p + q)\ip{0)j • nip(z) T 3 \7r a (p))\ + _ 



z+=0,z ± =0 ' 



The coordinate z lies on the light cone. At the quark- 
model scale the gluon GPD of the pion vanishes, 



7i (x, C, t) — (quark — model scale). 



(14) 



When the QCD evolution to higher scales is performed, 
a non- vanishing gluonic contribution H G (x,(,t) is gen- 
erated [Hj]. 

In the following analysis we use the symmetric notation 
for the GPDs, 



c 



2-C 



X 



x-(/2 
l-C/2' 



(15) 



where < £ < 1 and the support is — 1 < X < 1. One 
introduces the corresponding GPDs 



H I ~°' 1 (X,^,t) = H 



7=0,1 



t + X 2£ 



(16) 
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The reflection about X — yields the symmetry relations 



1=0/ 



H J ='(X,Z,t) = H^{-X^,t). 
At X > the GPDs are related to the PDF, q(X), 



(17) 



i? J=u (X,0,0) = iJ J=1 (X,Q,0) = q(X). 

The polynomiality conditions 0, Q follow from very 
basic field-theoretic assumptions, namely the Lorentz in- 
variance, time reversal, and hermiticity, which yields the 
form ([8|). For the moments one finds 

,1 3+1 

/ x 2j+l h i=o (x, e, o = 2 5^ A 

1=0 



2 ]+ 2,2,(i)f, (18) 



with j = 0, 1, ... . For the lowest moments one has 

JdXH I=1 {X^,t) = 2A w {t) = 2F v (t), (19) 
r -°(x^,t) ^ 2A 20 (t) + 2A 22 (t)^ 2 

= e 2 (t)-e 1 (i)£ 2 . (20) 



dXXH 1 
-1 



In the convention of Refs. 0, H[ the equivalent expansion 
is in powers of (2£) 2 rather than £ 2 , as in Eq. (fl~8)) . In our 
approach, used in the following Sections polynomiality is 
explicitly manifest from the use of the double distribu- 



tions [15[ . The equivalence of Eq. ([5]) and (fT5| is easily 
proven by contracting ([8]) with the null vectors n^ 1 . . . 
and subsequently applying the definitions (|11|) . We no- 
tice that for the isovector GPD only the even, and for the 
isoscalar GPD only the odd moments are non- vanishing. 
The gluon form factors are defined as the integrals 

f dx x 2 ^h 3 (x, e, t, q 2 ) = 2 J j2 A% +%2i (t) e. 

J - 1 i=0 

(21) 

Finally, we note that Eq. (fT8|) for j = expresses the 
electric charge conservation and the momentum sum rule 
operating in deep inelastic scattering. 



III. CHIRAL QUARK MODELS 

In chiral quark models at the leading- N c level the cal- 
culation of the form factors and GPDs proceeds according 
to the one-loop diagrams of Fig. [T] Extensive details of 
the quark-model evaluation are given in [ToT ]. In this pa- 
per we carry out calculations in two chiral quark models: 
SQM Q and NJL with the Pauli-Villars regularization in 




(a) 



p+q 



(b) 




p+q 




FIG. 1: (Color online) The direct (a), crossed (b), and contact 
(c) Feynman diagrams for the quark-model evaluation of the 
GPD of the pion. The contact contribution is responsible for 
the D-term. 



the twice-subtracted version of Refs. [4l|, |6J, |65(. Vari- 
ants of chiral quark models differ in the way of performing 
the necessary regularization of the quark loop diagrams. 

The spectral quark model 0] introduces the general- 
ized spectral density p(w) in the quark mass to, in the 
spirit of Ref . (66j , supplied with chiral symmetry, gauge 
invariance, and vector meson dominance. The one-quark- 
loop action of SQM has the form 



SQM 



= -iN, 



j dujp(uj)Ti-log(i$-ujU 5 
Jc 



(22) 



where p(u>) is the quark generalized spectral function, and 
U 5 = exp(«75T • <p/f), with cj) denoting the pion field in 
the nonlinear realization. The vector part of the spectral 
function, needed in the present analysis, has the vector 
meson-dominance form Ujl 



1 1 



1 



27r^(l-4w 2 /™ 2 ) 5/2 ' 



(23) 



exhibiting the pole at the origin and cuts starting at 
±m p /2, where m p ~ 770 MeV is the mass of the rho me- 
son. The complex contour C for the integration in 
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is given in Ref. 4]. SQM leads to conventional and suc- 
cessful phenomenology for both the pion d, [l5|, [62], [13] , 
the nucleon [68[ , and the photon parton distribution am- 
plitude 

We also study a more conventional chiral quark model, 
the NJL model with the Pauli-Villars regularization in 
the twice-subtracted version proposed in fiefs. (4f|, [64|, 
l65j . The one-quark- loop action of the model is 



Tnjl - -iN c Tr log (i@ - MU 5 ) 



(24) 



where M is the constituent quark mass. The Pauli- 
Villars regularization is introduced at the effective action 
level [4l|, IH, [65[ , with the practical advantage that gauge 
and relativistic symmetries as well as sum rules are man- 
ifestly fulfilled [70]]. For the observables considered in 
this paper the Pauli-Villars is implemented according to 
the prescription, where instances of M 2 in an observable 
O are replaced with M 2 + A 2 , and then the regularized 
observable is evaluated according to the prescription 



O(0) - C(A 2 



A 



2 dO(A 2 ) 
dA 2 ' 



(25) 



The Pauli-Villars regulator A is a free parameter of the 
model. In what follows we use 



M = 280 MeV, A = 871 MeV, 



(26) 



which yields / = 93.3 MeV for the pion decay constant 
[ill ] according the the formula 

Z 2 = -^[l°g(A 2 + M 2 )] rcg . (27) 



IV. ELECTROMAGNETIC AND 
GRAVITATIONAL FORM FACTORS 

The form factors may be calculated in two different 
ways. The first method uses the definition ([5]), which 
leads to the evaluation of one-loop diagrams with an ap- 
propriate vertex. For the electromagnetic form factor the 
vertex is Qj 1 *, where Q is the electric charge of the quark. 
For the gravitational form factor the vertex, correspond- 
ing to the energy-momentum tensor, has the form 

Q^(k + q, k) = i [(2k + qf^ + (2k + q)"<f] 

-l g ^ (2^ + 4- co), (28) 

with ijj denoting the quark mass. We illustrate the cal- 
culation in the Appendix. The other method uses the 
GPDs obtained earlier [l5[ and evaluates their moments 
([18]) . The results are the same, which serves as a consis- 
tency test of the algebra. As mentioned in Sec. HU the 
equivalence is proven by contracting with the null vector. 
For instance in the case of the gravitational form factor 
we consider n M 6 MI/ n„. Then 



(7r b (p + q ) |n M e^(0K Kb)) = 
<5 afc i[c 2 e 1 (q 2 ) + (2-o 2 e 2 ( g 2 )] 



(29) 



and the vertex becomes 

n li & a '(k + q,k)n v = (x-C/2)>yn. (30) 

We recognize the same vertex as in the evaluation of the 
GPD's multiplied by (x — C,/2). Upon passing to the 
symmetric notation Eq. (|20[) follows. 

Before showing the explicit results both for SQM and 
NJL models in the specific realizations described in the 
previous section, we note some general results. Actually, 
in the considered quark models and in the chiral limit we 
have the following identity relating the gravitational and 
electromagnetic form factor, 



j f [tem = F v (t), 



(i = l,2) 



(31) 



from which the identity between the two gravitational 
form factors 9i(i) = 9 2 (i) = <d(t) follows. This re- 
markable quark model relations are proven explicitly in 
the Appendix. The essential ingredient of the proof is 
the existence of the spectral representation in both con- 
sidered models. One consequence of relation (f3"Tj) is the 
expected consistency of normalizations at t = Q for the 
charge and mass 6^(0) = i*V(0), displaying the tight 
connection between the gauge and Poincare invariances. 
Furthermore, expanding in small t and using the fact 



that F(t) = F(0)[1 - (r 2 )t/6 + . . 
the mean squared radius, one gets 



2(r 2 ) e 



(r 2 )v, 



with (r 2 ) denoting 



(32) 



which means that in the considered models and in the 
chiral limit the gravitational mean squared radius is half 
the electromagnetic one. 

Turning now to the specific realization, by construc- 
tion, the pion electromagnetic form factor in SQM has 
the monopolc form 



F* QM (t) 



m 2 p -t' 



(33) 



A straightforward evaluation for the gravitational form 
factor yields 

Qf^(t) = e^ QM (t) = ^ log = Q SQM (t). 

(34) 

These specific form factors fulfill trivially the general re- 
lations Eq. J3ID and Eq. (j52")l . 

In the NJL model the pion electromagnetic form factor 
is equal to 



pN.TL 



(t) 



Nc9l 
8tt 2 



2 s log 



s+V-t 



with = Mj f denoting the quark-pion coupling con- 
stant. We have introduced the short-hand notation s = 



6 



Aio 




FIG. 2: (Color online) The electromagnetic form factor (top) 
and quark part of the gravitational form factor (bottom) in 
SQM (solid line) and NJL model (dashed line) compared to 
the lattice data from Ref. The band around the SQM re- 
sults corresponds to the uncertainty in the quark momentum 
fraction R and the m p parameter. 



^4 (M 2 + A 2 ) -t. The condition limt^-oo F NJL (t) = 
is satisfied due to Eq. (|27p . For the gravitational form 
factor we obtain 

9 NJL W = ^NJL^ _ M x (35) 

where hi2(z) is the polylogarithm function. 

Of course, the explicit expressions ()33I34|) and (|35I35|) 
comply to the general relation (|3"Tj) . 

Similar expressions, of growing complication, may be 
obtained for higher-order form factors in both SQM and 
NJL models. 



V. COMPARISON TO THE LATTICE DATA 

We are now ready to compare the results of the pre- 
vious Section to the recent lattice data @, These 
data correspond to the scale of Q — 2 GeV and give 
the quark part of the gravitational form factor and the 



-tAio [GeV 2 ] 



0.8 r 




-tA 20 [GeV 2 ] 



0.5 




FIG. 3: (Color online) Same as Fig. [2] for the form factors 
multiplied with —t. In addition, we include the TJLAB data 
fFH . [72 . [73| ] (darker and larger circles and squares) and the 
Cornell data [7|| (diamonds). 



electromagnetic form factor. On the other hand, the 
quark model calculation corresponds to a very low scale, 
Qo ~ 320 MeV and in general the QCD evolution is 
needed to compare the model predictions to the data at 
a different scale Q. 

A detailed discussion of the evolution issue is presented 
in Ref. [15] . The low energy Chiral Quark Models contain 
quarks as the basic and explicit degrees of freedom and 
the considered twist-2 observables correspond to model- 
ing QCD at a low renormalization point. Obviously, by 
the energy-momentum conservation in these models the 
quarks as the only degrees of freedom carry 100 % of the 
total momentum in a hadron. On the other hand the mo- 
mentum fraction carried by the valence quarks in the pion 
at the scale Q 2 = 4GeV 2 is about 40%. The QCD evolu- 
tion (LO, NLO, N2LO) tells us that this number grows 
as the scale is evolved to lower values. The quark model 
reference scale Qo is determined as the scale where the 
evolved QCD value yields exactly the inescapable 100% 
of the quark model. It would of course be highly desir- 
able to include, e.g., explicit gluonic degrees of freedom 
as this would allow to stop the perturbative evolution at 
a higher scale. Nonetheless, the LO and NLO evolutions 
for the PDF's are sufficiently close to each other [l8[ as 
to provide some confidence on the kind of calculations 
carried out in our work. 



7 



The issue of the QCD evolution of GPDs and the gen- 
eralized form factors is addressed in detail in Sect. IVII 
However, the electromagnetic and gravitational form fac- 
tors do not evolve with the scale. What changes is the 
ratio R of the total momentum fraction carried by the 
quarks (valence and sea) 



R = 



) q (Q) 



a(Q) 



7i 0) /(2/3o) 



(36) 



where the anomalous dimension is given by 7x / (2/3o) 
32/81 for N F = N c = 3 and 



a(Q 2 ) 



-in 



/3 log(Q 2 /A 2 



QCD ) 



p = ^N c - 2 -N f , 



(37) 
(38) 



where we take A QC d = 226 MeV and N c = N f = 3. 
At the scale Qo we have R = 1, which then gradually 
decreases with the increasing scale. The quark part of 
the gravitational form factor is Q q (t) = RQ(t), the gluon 
part is e G (t) = (l-i?)9(t), such that, of course, @ q (t) + 

e G (t) = e(t). 

The value of R depends on the evolution ratio 
a (Q 2 ) I a (Qo) i which is not precisely known on the lattice. 
For that reason we shall treat R as a free parameter when 
fitting the model results to the data for the gravitational 
form factor. In SQM the other parameter is the value 
of the rho meson mass, m p . For SQM we fit jointly the 
electromagnetic and the quark part of the gravitational 
form factor. The \ 2 method yields 



R = 0.28 ±0.02, 



to p = 0.75 ± 0.05 GeV, (39) 



with x 2 /DOF = 1.8. The result of the fit is displayed 
in Fig. [2] with the solid line. The band corresponds to 
the uncertainties in the values of parameters in Eq. (|39p . 
We note an overall very good agreement for both form 
factors. Note a significantly slower fall-off for the gravita- 
tional form factor compared to the electromagnetic one. 
The optimum value of m p agrees within the error bars, 
which are substantial, with the physical mass of the rho 
meson. 

For the case of the NJL model (dashed line in Fig. [2|) 
the agreement is not as good as in the SQM model and it 
is not possible to improve it by changing the parameters 
M and A. For that reason we have not carried out the 
X 2 fit in this case. The problems of the NJL model in 
reproducing the electromagnetic form factor, where the 
corresponding rms radius turns out to be too small, are 
well known, see e.g. the discussion in Ref. [ZD ]. Similar 
discrepancy can be noted for the case of the gravitational 
form factor shown in the bottom panel of Fig. [51 

Now we come back in a greater detail to the issue of 
the quark to gluon momentum ratio R. It evolves with 
the scale from the value R = 1 at the quark model scale 
to R -> at Q 2 -> oo. We use the standard LO DGLAP 



evolution for the single channel. We know the final scale, 
Q = 2 GeV, but the quark model scale is a priori un- 
known. We thus adjust Qo m order to reach the value of 
R from Eq. (f3"9"|) at the known scale Q. The result is 



Qo = 0.31 ± 0.03 GeV 



(40) 



with Aqcd — 226 MeV. This value agrees with the ear- 
lier independent determinations of the quark- model scale, 
see Ref. [15j for a detailed discussion on these issues. 
We note that at such low scales as in Eq. (|4"0"]) the per- 
turbative expansion parameter in the DGLAP equations 
is large, a(Q^)/(2ir) — 0.34, which makes the evolution 
very fast for the scales in the vicinity of Qo- 

One should note that the large value of a(Qo) calls for 
the use of improved QCD evolution at low scales. The 
presented analysis could be extended in several ways, for 
instance incorporating the next-to-leading (NLO) effects, 
or by modifying the dependence of a on Q 2 at low scales 
incorporating the "infra-red protection" [7a |. The NLO 
corrections to the meson PDF's in chiral quark models 
were studied in Ref. [llj], where, somewhat surprisingly, 
it was found that these effects are small. Modifications in 
the evolution of ct(Q 2 ) co uld be incorporated along the 
lines of Refs. [7y, u3, lZ8| where they were used for the 
pion form factor, but this calculation is outside of the 
scope of the present paper. 

We stress that the chiral quark model explanation 
of the lattice data used in this work as well as the 
physical and lattice data explored in Ref. [l5T ] requires 
"strong" evolution, with a large value of the evolution 
ratio a(Qo)/a(Q 2 ). Our simple-minded analysis based 
on the LO evolution may probably be viewed as an ap- 
proximation to a more elaborate scheme, nevertheless it 
is quite remarkable that various observables (PDF, PDA, 
their moments) lead to the same evolution ratio, giving 
in the LO approximation compatible quark-model scales 
flU . Moreover, the electromagnetic and the gravitational 
form factors are independent of the evolution, hence the 
issue does not arise for these observables. 

In Fig. [3] we show the form factor multiplied by — t, 
which is a popular way to present the results at large 
Euclidean momenta. In the case of the vector (electro- 
magnetic) form factor we also display the experimental 
TJLAB data [lE Iz3, III and the earlier Cornell data 
[zl] and the approximately constant value for —tFy{t) 
is clearly seen. In the gravitational case, however, the 
lattice data 0, 0| show an increasing trend which is well 
mimicked by our SQM form factor, —tO(t) ~ log(— t). 

In the large- N r . limit in the single resonance approxi- 
mation (SRA) [6l|, |79j one would have a monopole form 
factor 



e BKA (t) 



AM] 
MfTt 



(41) 



which agrees with our quark model at low t if Mf = 
y/2M v ~ HOOMeV in agreement with the result of 
Eq. (f3"2"l) . The difference between the SRA form factor 
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and our Eq. (f34|) is less than 10% for momentum val- 
ues up to t = — lGeV 2 . The SRA monopole fits well 
the data extrapolated to the physical pion mass, yielding 
A = 0.261(5) and M f = 1320(60)MeV 

To conclude this Section, several general remarks are 
in order. One has to bear in mind that the full-QCD lat- 
tice calculations are linearly extrapolated to the physical 
pion mass. Our quark-model analysis incorporates only 
the leading- N c contributions, while the full-QCD simula- 
tions include all orders in that expansion. Nevertheless, 
despite these caveats, we note a quite remarkable agree- 
ment, in particular for the case of SQM. This is notewor- 
thy, as our study is a genuine dynamical field-theoretical 
calculation, and the parameters providing the optimal fit 
are highly compatible with calculations of other processes 
within the same model and scheme. 

Admittedly, at very large values of —t perturbative 
QCD results for the form factors should be reached. At 
very low values of t chiral corrections are important. 
Since the data we use are at intermediate values of Eu- 
clidean momenta, we may neglect both above-mentioned 
effects and use the chiral quark models to explain the 
data. 

We stress that our calculation conforms to the low- 
energy theorem 0i(O) — 6a(0) = 0(m^) [6l[ which is 
dictated by chiral symmetry. The quark model predicts, 
in addition, Oi(i) = Q 2 (t), and the purely multiplicative 
character of the QCD evolution yields in our conventions 



A 22 (t) = -A 20 (t). 



(42) 



Probably due to insufficient statistics, this formula is not 
quite seen in the data of Ref. Q (note that with conven- 
tions adopted in that reference one should have instead 
^4-22 (t) = — jA 20 (t)). For that reason we have not used 
the data for A 22 in our numerical analysis. 



VI. HIGHER-ORDER FORM FACTORS 

Analogous calculations as in the previous section can 
be performed for the higher-order generalized form fac- 
tors. Here we only give the results for the case of SQM, 
as the results in NJL are qualitatively similar, while also 
SQM works better for those quantities which can be con- 
fronted to the data. In Fig. Ufa) we show A 3y2 i and A^ 2i 
obtained at the quark model scale. In the chiral quark 
models in the chiral limit one has at t = very simple 
expressions [H, ED, HI 



H J=i (X,£,0) 



l-X A 



(43) 



H 



1=0 



e((i - x)(x -0)-e((x + -x)), 



the following relations hold [5j 




2 j+ i forz = 
otherwise 

for i = 

2j+2 fori = j + l 
otherwise 



(quark — model scale) (44) 

These relations can be seen in Fig.fJJa). The form factors 
tend to zero very slowly at large — t. 

Another property follows from the fact that in the con- 
sidered model H I=0 (X, l,t) = for any value of t. Then 
Eq. HID yields [H 



j'+i 

i=0 



i(t) = o. 



(45) 



This feature can be seen in the lower panel of Fig. |U 

While the lowest order form factors, i.e. the electro- 
magnetic and the gravitational, do not evolve with the 
running scale Q, the higher-order form factors change 
when we pass from the reference quark model scale Qo 
to Q. Our analysis proceeds as follows: We carry out 
the leading-order DGLAP-ERBL evolution for the quark- 
model GPDs from the quark-model scale 



Qo = 0.31 GeV 



(46) 



to the scale of the lattice calculation, Q — 2 GeV. Note 
that the GPDs at the quark-model scale do not exhibit 
factorization in the i-variable. In fact, this is the rea- 
son for the non-trivial change of the higher-order form 
factors with the scale. After the evolution we take the 
moments (fT5j) at several values of £ and then disentangle 
the generalized form factors via solving a set of linear 
equations. For the LO ERBL DGLAP evolution we use 
the method and the numerical program of Ref. [80( . The 
result of our calculations is presented in Fig. 0Jb). We 
notice a dramatic change, both in the value at t = and 
in shape of the form factors compared to the behavior of 
Fig. [Ua) at the quark-model scale. 

Formally, as Q 2 — > oo the GPDs approach their asymp- 
totic forms contained entirely in the ERBL region |A| < 
£. Explicitly, we have in this limit [15J 



H 



i=i 



3 

2i 



1 



H 



1=0 



XH 



G 



d-e) 

(Q 2 



A 2 

e 

15 



F v (t) 



N 



I 



X 



4£ 2 4C F + N f £ 
15 C F 



4£ 4C F 



N 



1 



e 



(47) 

e(t) 



— ) e(*)> 



where 9(x) is the Heavyside step function. It follows 
from the definition (fT5|) that at the quark-model scale 



where C F = (N 2 - l)/(2iV c ) and N f = 3. The propor- 
tionality factors follow from the normalization at the ini- 
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FIG. 4: (Color online) Generalized form factors ^.3,21 and 
Ai^i of the pion in SQM at the quark-model scale Qo (a), at 
the lattice scale Q — 2 GeV (b), and the gluon form factors 
Af 2i at Q = 2 GeV(c). 



tial quark-model scale Qo, as the charge- and momentum- 
conservation sum rules are invariants of the evolution, 



J dXH I = 1 (X > £,t,Q 2 ) = 2F v (t), (48) 
J ^ dX [XH I=0 (X, £, t, Q 2 ) + XH g (X, £, t, Q 2 )] 

in accordance to Eq. (|19I20[) . 



Evaluation of moments in Eq. (|47[) yields immediately 
r 3 F v (t) fori = j 



A^2j+l,2i(t) 
A2j+2,2i(t) 



4i(i+2)+3 J 

otherwise 

4C F +N f 2[4j(j+4) + 15]®(*) forl =J 

-^■2i+2,2j(*) for i = j + 1 

otherwise 



4' 



(t) 



4C 



1 



^V> 2j + l" 



-.4 



2j+2,2i 



(Q 2 



00 



(49) 



We note a qualitative difference of the asymptotic form 
factors compared to the form factors at the quark model 
scale shown in Fig. [4] For the isovector case (n = 2j + 1) 
only the highest form factor, with i = j, does not vanish, 
while for the isoscalar case (n = 2j + 2) only the two high- 
est moments, with i = j + 1 and i = j are non-zero. The 
remaining moments tend to zero. This result is a prompt 
conclusion from the asymptotic forms (|4"T)) . Note that 
in contrast to this behavior, at the quark-model scale 
all generalized form factors are non-zero. As mentioned 
previously, the form factors A w , A20, and A22 are invari- 
ants of the evolution. The asymptotic gluon form factors 
in Eq. (|4"9"|) are related to the isoscalar quark form fac- 
tors in a simple manner. Asymptotically, all generalized 
form factors become proportional to Fy(t) or <d(t) in the 
isovector and isoscalar channels, respectively. 

Finally, we compare our values of the higher-order form 
factors at t = to the lattice data provided in Sec. 7 of 
Ref. 3. With the notation for the moments at t = 0, 

(x n ) = A n+1>0 (0), (50) 

one finds at the lattice scale of Q = 2 GeV 

(x) = 0.271 ± 0.016, (51) 
(a; 2 ) = 0.128 ±0.018, 
(a; 3 ) = 0.074 ±0.027. 

(lattice) 

while in both chiral quark models we obtain after the LO 
DGLAP evolution to the lattice scale 

(x) = 0.28 ± 0.02, (52) 
(x 2 ) = 0.10 ±0.02, 
(x 3 ) = 0.06 ±0.01, 

(chiral quark models) 

where the error bars come from the uncertainty of the 
scale Qo in Eq. (|4"0)) . The two sets of numbers overlap 
within the error bars. This result is quite remarkable, as 
it shows that the hierarchy of the form factors at t = 
found in full-QCD lattice calculations is properly repro- 
duced in chiral quark models. The above form factors 
are simply the moments of the PDF of the pion. We 
recall that the PDF itself in the chiral quark models re- 
produces very well the experimental [52 1 and transverse 
lattice data of Ref. [53| . see Figs. 8 and 9 of Ref. 
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VII. CONCLUSIONS 

Here are our main points: 

1. The gravitational form factor obtained from the 
spectral quark model in the chiral limit agrees very 
well with the lattice data of Refs. @, ||. We have 
performed a global fit to the electromagnetic and 
the quark part of the gravitational form factors and 
obtained very reasonable values for the evolution 
ratio and the vector-meson mass. The longer range 
of the gravitational form factor follows from a dif- 
ferent analytic expression compared to the electro- 
magnetic form factor. 

2. The NJL model does not provide such an excel- 
lent agreement, although the qualitative features 
are very similar to SQM. 

3. We provide analytic expressions for the lowest- 
order form factors in both considered models. For 
the considered models in the chiral limit we find an 
explicit relation between the gravitational and vec- 
tor form factors. In particular, the relation shows 
that in our case both gravitational form factors are 
equal, and that the mean squared electromagnetic 
radius is twice the gravitational one. 

4. The electromagnetic and gravitational form factors 
do not evolve with the scale, but the higher-order 
generalized form factors do. We have performed 
the leading-order ERBL DGLAP QCD evolution 
of the pion GPDs and obtained via moments the 
generalized form factors at the scale of the lattice 
measurements. 

5. The generalized form factors at t = found in full- 
QCD lattice simulations are reproduced in chiral 
quark models within the error bars corresponding 
to statistical and model uncertainties. 

6. Our predictions can be further tested with future 
lattice simulations for higher-order form factors. 
The behavior is non-trivial, with form factors hav- 
ing different signs, magnitude, and asymptotic fall- 
off. 

7. Lattice simulations for the gluon form factors would 
provide a very useful independent information, 
which could be used to verify the model predic- 
tions. 
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APPENDIX A: RELATION BETWEEN 
GRAVITATIONAL AND ELECTROMAGNETIC 
FORM FACTOR 

In this Appendix we show an explicit calculation of 
form factors in SQM and prove as a by-product the result 
of Eq. (|3Tj) . The proof also applies to the NJL model with 
the PV regularization, as that model can also be written 
in terms of the spectral representation (see below). 

The SQM is defined by the generalized Lehman repre- 
sentation of the quark propagator 



S{j>) = I duj 
Jc 



p(uj) 



OJ 



(Al) 



The method exploited in Ref. [|| was mainly based in 
writing down Ward identities for the corresponding ver- 
tex functions at the quark level, and then closing the 
quark line to make one loop calculations. As discussed 
in later works (see e.g. [H, [68|) this is fully equiva- 
lent to proceeding directly through the effective action 
which we sketch now. Using the action of Eq. (f2"2"|) one 
can compute the energy momentum tensor as a func- 
tional derivative of the action with respect to an exter- 
nal space-time-dependent metric, g^(x), around the flat 
space-time metric r/^ (we take the signature (H ) ) 



2 U [X) Sg^(x) 



(A2) 



i^ f dujp(oj)(x\{o^,^-ojU 5 ) 1 }\x), 



where 



= - {^d v + -y v d») - of (i# - w) , (A3) 

and U 5 — exp(i75T • 4>/ f) with </> the pion field in the 
nonlinear realization. For the calculation of the gravi- 
tational form factor, we expand up to second order in 
the pion field corresponding to evaluating the diagrams 
of Fig. [T] In the cartesian isospin basis one has 



(^le^K) = -n c / du P (u) 



d A k 



rTr 



- - \ f 

l iS a bUJ 



— j) — OJ 

i 



(2tt)4 

/ 

crossed| , 



where the quark gravitational vertex is given by 



(A4) 



1 r 

4 
1 



{2k ll + q^-fv + (2k u + q v )^ p 



g^fi + i-uj) 



(A5) 



After computing the traces and using the Feynman trick 
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in the integrals, in the chiral limit the result is 



tor reads 



e^ 2 



4tt 2 / 2 



du> p(cj)u' 



dx 



x(l~x)q 2 



- ui 2 logw 2 + x(l - x)q 2 
+ {u 2 - x{l - x)q 2 ) log(w 2 - x{\ - x)q 2 ) 

e 2 ( 9 2 ) - e^ 2 ). 



(A6) 



In the low-momentum limit we may use Eq. © as de- 
duced from xPT in the presence of gravity [6l| to get 



in 



Nr. 



L\2 — — 

L13 = — 



1927T 2 

Nr. 



96tt 2 ' 

{An) 2 12B 2AM 2 ' 



(A7) 
(A8) 
(A9) 



in agreement with the derivative expansion of SQM [62J. 
Proceeding similarly with the vector form factor we get 



(^ a | J£'V) = -N c / dwp{u) 



Tr 



d A k 
(2tt)* 



n 

7 "2 



For on-shell massless pions the electromagnetic form fac- 



(A10) 



iM<z 2 ) = - 



du>p((jj)uj 2 



Air 2 f 2 
l 

x / da; log [to 2 - x(l - x)q 2 ] 



(All) 



Charge normalization, Fy(0) — 1, and energy- 
momentum normalization 02(0) = 1 implies 



i = f v (o) = e 2 (o) 



4tt 2 / 2 
N c m 2 p 
24tt 2 / 2 



dujp(u>)cj 2 log w 2 



(A12) 



where in the second line the vector meson realization, 
Eq. (f2"5| . has been used. The value agrees with the value 
of the pion weak decay constant obtained from the cor- 
responding axial matrix element 4]. 

With the representations of Eq. (|A6|) and Eq. (|A1 1[) 
the result in Eq. (|31[) can be readily derived, without any 
reference to the specific realization given by Eq. (|23|) . 

The above proof also holds for the NJL model in 
the PV regularization. This class of models can be 
cast explicitly in the spectral-representation form, using 
p(uj) = S(uj-M) + J2i CiS(u- v/M^+Af), where c t and 
Aj are the PV constants. Then all above algebraic steps 
carry over and the result (|3Tj) holds for the NJL model 
with PV regularization as well. 
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